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Abstract 

Recent measurements of the branching fraction for 
D s — > Iv disagree with the Standard Model by around 
2ct. In this case the key aspect of the Standard Model is 
the calculation of the decay constant. fr> s , with lattice 
QCD. This talk surveys the experimental measurements, 
and explains how the lattice QCD calculations are done. 
Should the discrepancy strengthen again (it was earlier 
3.8er), it would be a signal of new physics. Models that 
could explain such an effect are also discussed. 

1. Introduction 

The decay constant of a meson parametrizes the decay 
of the meson to leptons. For a pseudoscalar like the 
charmed strange meson D s , it is defined by the hadronic 
matrix elements 

(0|s7 M 7 5 c|£> s ) = if Ds P^ (1) 
(m c + m s )(0\sj 5 c\D a ) = -m 2 D Jr> 3 . (2) 

The Ward identity of the partially conserved axial cur- 
rent (PCAC) ensures that the two definitions are iden- 
tical. These matrix elements are directly computable 
in QCD, via numerical simulations of lattice gauge the- 
ory. These calculations are useful: the ratio f^/fx, 
for example, is used to determine the Cabibbo angle, 
tanflc c< [B(K lv)/B(n lv)Y' 2 U/f K [I]. 

In the Standard Model, the expression for the branch- 
ing ratio is 

B(D S -+ tu) = (l ~ £f) 2 l<fcO*| a , 

(3) 

with analogous expressions for leptonic decays of other 
pseudoscalar mesons. Here Gf is the Fermi constant, 
measured in muon decay, and V is the CKM matrix. 
Beyond the Standard Model, one replaces 

2 

G F V*m e h-> G F V*m e + G A m e + Gp m ° B , (4) 

m c + m s 

where G A and Gp are related to couplings and masses of 

new interactions in a way analogous to Gf = g 2 / V^M^- 
Amplitudes that proceed through an axial- vector current 
(Gf and Ga) are helicity suppressed, but amplitudes 
that proceed through a pseudoscalar current (Gp) are 
not. In practice, D s — ► ^lv has a helicity-suppression 
factor m 2 ll /m 2 Ds = 2.8 x 1CP 3 . On the other hand, D s — > 
tv is not helicity suppressed but, instead, phase-space 
suppressed: (1 — m 2 /m 2 Da ) 2 = 3.4 x 10~ 2 . 

The decay constant of the D s was expected to be an 
excellent test of lattice QCD, for several reasons [2]. The 
matrix elements in Eqs. (UJ and © are gold-plated, in 
the sense of Ref. [3], namely, only one hadron enters, 
and the chiral extrapolation is controlled. Experimen- 
tal measurements of |V r cs |/D s , via Eq. can be com- 
bined with the determination of \V CS \ from CKM unitar- 
ity. The idea that new physics could compete is usually 
discounted, because the decay is Cabibbo-favored and 



proceeds at the tree level of the weak interactions. Fi- 
nally, the precision of experiments has lagged that of cal- 
culations, so the analysis of the numerical lattice-QCD 
data is carried out with a relatively blind eye. 

The first round of testing seemed to go well. In June 
2005, the first lattice-QCD calculation with 2+1 flavors 
of sea quarks appeared in a joint work of the Fermilab 
Lattice and MILC Collaborations: 

f Ds = 249 ± 3 ± 16 MeV Fermilab/MILC g], (5) 

where the first error is statistical and the second system- 
atic. This prediction was followed a year later by a com- 
parably accurate measurement of B(D S — > jiv) / B(D a — > 
07r) that, when combined with an independent measure- 
ment of B(D S — ► 07r), yielded 

f Ds = 283 ± 17 ± 7 ± 14 MeV [iv/fa BaBar [5], (6) 

which agrees with Eq. ([5]) at 1.2er. Both results are dis- 
cussed in more detail below. 

But then something unexpected happened. During 
2007 the CLEO and Belle Collaborations both pub- 
lished absolute measurements of B(D S — > (M/) [SJ [7], 
and CLEO also published absolute measurements of 
B(D S — > tv) [HI IB]- Transcribed via Eq. as mea- 
surements of the decay constant, these are 

f D , = 264 ± 15 ± 7 MeV fxv CLEO [S], (7) 

f Dg = 275 ± 16 ± 12 MeV fjtu Belle [7], (8) 

Jd s = 310 ± 25 ± 8 MeV tv CLEO [5], (9) 

f Ds = 273 ± 16 ± 8 MeV tv CLEO B]. (10) 

Taking Eqs. (|6"))- (fl"t)|) at face value, the weighted average 
(combining all errors in quadrature) is [9] 

fo, = 277 ± 9 MeV, (11) 

which is 1.5cr higher than the value in Eq. Mean- 
while, and more dramatically, the HPQCD Collabora- 
tion published a lattice-QCD calculation with an error 
significantly smaller than Fermilab/MILC's: 

f Ds = 241 ± 3 MeV HPQCD [TO]. (12) 

Sect ion [3] explains why t he e rror i s so much smaller. The 
difference between Eqs. (fTTj) and (fl"2"|) is 3.8ct. (Omitting 
BaBar, as in Ref. [11], the discrepancy becomes 3.4(7.) 
It is imp ortan t to bear in mind — and easy to see from 
Eqs. (fTTj) and (fT2|) — that the yardstick for a is the exper- 
imental statistical error. The 2008 edition of the Review 
of Particle Physics [TTJ noted that the discrepancy could 
be a sign of physics beyond the Standard Model. Can- 
didate models are discussed below. 

The rest of this paper brings this saga up to date. Sec- 
tion [23 gives a brief survey of the experiments, including 
higher-statistics measurements from CLEO, and a (dif- 
ferent) extraction of Jd 3 from BaBar's measurement of 
B(D S — > fiv)/B(D s — > <j>ir) by the Heavy Flavor Averag- 
ing Group (HFAG), which have brought the discrepancy 
down to 2.3(7. Section [5T] discusses recent developments 
in lattice-QCD calculations. New physics explanations 
are in Sec. 23 The main issues are summarized in Scc.[5~1 
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2. Measurements 

Observations of D s — ► £v date back to 1993 in fixed- 
target experiments and e + e~ collisions [12j [13j HH [T5l 
[TBI 23 [HI [19j EDJ. These early measurements are omit- 
ted from (the current) Particle Data Group (PDG) and 
HFAG averages and are also omitted from this discus- 
sion. In 2006, the PD G |2"T ] i ncluded a correlated av- 
erage of Refs. [IDJ OH [T8 p T9 l [2D] , which increase the 
discrepancy, as discussed below, by 0.3-0.4er. 



2.1. CLEO pv and tv 

(*) 

CLEO produces D S D S pairs in e + e~ collisions just 
above threshold, as in a 1994 observation by the BES 
Collaboration [Ti] . The multiplicity is low, so the whole 
event can be reconstructed, and the neutrino is "de- 
tected" by requiring the missing mass-squared to be con- 
sistent with 0. Radiative events with photon energy 
greater than 300 MeV are rejected. Although this cut is 
imposed for other reasons, it usefully removes radiative 
events without helicity suppression. 

In D s — > tv, the r decays in the detector, and the de- 
tails of the analyses depend on the r-decay mode. CLEO 
first observed events in which r — > %v as a background 
to the D s — > fxv analysis, but then turned these events 
into a measurement. A separate analysis chain counts 
D s — > tv events in which r — > evv. With 2 or 3 neu- 
trinos in the final state, the constraint on missing mass- 
squared is no longer pertinent. These analyses also reject 
events with photons, but this is a matter of r detection. 
In the D s rest frame, the r acquires only 9.3 MeV of 
kinetic energy, so radiative events are not an issue. 

In January 2009, CLEO published analyses with their 
full data-set, reporting 

fo, = 257.3 ± 10.3 ± 3.9 MeV pv [22], (13) 
fn B =278.7 ±17.1 ±3.8 MeV tv, t -> nv [22J, (14) 
f D , = 252.5 ± 11.1 ± 5.2 MeV tv, t -> evv [23], (15) 

which supersede Eq. ([7]), ^ and ([TP]) , respectively. After 
Physics in Collision 2009, CLEO made public an analysis 
of a third D s — > tv decay chain, r — > pv, yielding 

f D , = 257.8 ± 13.3 ± 5.2 MeV tv, t -> pv [24]. (16) 

A novelty of this analysis is that it disentangles a mesa- 
shaped signal distribution from a peaking background. 



2.2. BaBar and Belle \iv 

BaBar and Belle, following a strategy devised by 
CLEO [131 HZ], collect a D s sample from continuum 
events under the T(45) by observing the decay D* — ► 
D s j. BaBar then counts the relative number of D s j — > 
/iv"/ and _D S 7 — ♦ 07T7 events, yielding a measurement of 
B(D S — ► pv)/B(D s — > 4>tt). A separate measurement 
of B(D S — > (f>ir) is needed to extract fn, via Eq. ([3]), 
and B aBa r used an average of two of its own measure- 
ments [25]. Belle improves on the D* — > D s ^f technique 
by devising a Monte Carlo analysis to guide full recon- 
struction of the event. In this way they obtain an abso- 
lute measurement of B(D S — > pv). 

Measurements of B(D S — ► pv)/B(D s tfiir) are sub- 
ject to some ambiguity. The <fi decays to KK, but other 
processes, such as D s — > /o7r — ► KKir also contribute. 
The two contributions arc not completely separable, be- 
cause the amplitudes interfere [2D] . 



I?. 5. CKM; Radiative Corrections 

To extract /d s from the measurements of the branch- 
ing ratio, one needs a value of the CKM matrix ele- 
ment \V c d\- In practice, it has been determined from 
CKM unitarity, either using a global fit or simply set- 
ting 1 14a | = (R makes an insignificant difference.) 
With four or more generations, this assumption incor- 
rect, but 4- (or more) generation CKM unitarity still 
requires \V CS \ < 1. Therefore, an incorrect assumption 
about | V cs | cannot explain why the "measured" value of 
Jd s is too high. 

Leptonic decays are, of course, subject to radiative 
corrections. A class of virtual processes are of special 
interest here, namely D s — > D*7 — > pvj, where D* is 
a vector or axial-vector meson. The decay D* s — * pv is 
not subject to helicity suppression, so the absence in the 
rate of a factor (m^/mrj,.) 2 could compensate for the 
presence of the factor a f» 1/137. The radiative rate is 
significant for energetic photons [27] [28]. With CLEO's 
cut rejecting radiative events with E 1 > 300 MeV, how- 
ever, Eq. (12) of Ref. [27] shows that these events add 
only around 1% to the rate and, thus, cannot be an ex- 
planation of the discrepancy. 



2.1 HFAG 

The experimental collaborations' differing treatments 
of 1 14s | and of radiative corrections are not yet signifi- 
gant, so Refs. [9] [IT] simply average quoted values of 
fr>, ■ Eventually, however, a uniform treatment will be 
necessary, so, with this i n mind , the Heavy Flavor Av- 
eraging Group (HFAG) [2D] [3D] has undertaken to av- 
erage the model-independent quantities B(D S — > pv), 
B(D S -> tv), and B(D S -> pv)/B(D s -> <jnr). The av- 
eraging is straightforward. When turning to the extrac- 
tion of /d s , however, HFAG noticed an important issue 
with BaBar's determination of fo s ■ The definition of the 
4> resonance in Ref. [S] is a window of K + K~ invariant 
mass M K + K -, such that \M K + K - — < 5.5 MeV. 
The normalizing measurements, on the other hand, used 
\M K + K - - m,p\ < 15 MeV [25]. From the M K+K - dis- 
tribution in Ref. [2D], it is clear that the difference is 
important. Fortunately, CLEO [2DJ reports B(D S — ► 
A' + A'~7r) as a function of M K + K -, so HFAG combines 
B(D S -> K+K-tt) with \M K+K - -ro^l < 5 MeV, 
B((f> -> K+K-), and BaBar's B(D S -> pv)/B(D s -> 
4>n) to arrive at B(D S — > pv). Interpreting this branch- 
ing fraction as Jd b yields 

f D , = 237.3 ± 16.7 ± 1.7 MeV pu/cjm HFAG [3D], (17) 

which we shall use to supersede Eq. ([6j. It is 16% or 2.9cr 
lower (using the normalization and systematics for a). 



2. 5. Synopsis 

In summary, the measurements of the branching frac- 
tion B(D S — > tv) are relatively straightforward counting 
experiments. They can be contrasted with, say, searches 
for the Higgs boson at hadron colliders [3T], in which a 
careful and subtle modeling of the QCD background is 
essential. Here the background is small and-or measur- 
able; the events are clean, or even pristine. As a result, 
the dominant experimental error is statistical. It is, of 
course, possible that more experiments have fluctuated 
up than down. 

With the new results, including Eqs. (fTDJ and (TIT]) , 
the experimental average is now (I find) 



f Ds = 257.8 ±5.9 MeV, 



(18) 
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or 1.7a lower than Eq. (jlip . which is a combination of 
1.3cr from CLEO's new measurements and f.lcr from 
HFAG's revision of BaBar's measurement. (My aver- 
age without Eq. ([16]) is 257.8 ± 6.4 MeV, which is close 
to HFAG's more rigorous average of the same inputs, 
256.9 ±6.8 MeV [30].) 

3. Lattice QCD 

Lattice QCD has made great strides in the past several 
years [3] [35], compared to, say, the status at Physics in 
Collision 2002 [33] . The key development has been the 
inclusion of sea quarks, first with n/ = 2 and, then with 
?i/ = 2+1. The latter notation means that one sea quark 
has a mass nearly equal to that of the strange quark, and 
the other two vary over a range 0.1m s < m q < 0.5m s , 
such that chiral perturbation theory can be used to reach 
the up- and down-quark masses. 

That said, there have been only two calculations o f fp 
and Jd b with rif = 2 + 1 flavors of sea quarks [4] [10] , 
one of which dominates the average. Moreover, both use 
the same ensembles of lattice gauge fields [321 [34] , which 
have been generated using "rooted staggered fermions" 
for the sea qua rks. The rooting procedure leads to some 
difficulties [35] th at are expected to go away in the con- 
tinuum limit [35]. 

The reason for the rooting is that lattice fermion 
fields corresp ond to more than one species in the con- 
tinuum limit [37] , With staggered fermions there are 
four species [38] 7 Sea quarks are represented by a de- 
terminant of the (lattice) Dirac operator, so to reduce 4 
species to 1, one can make the Ansatz [35] 



Table 1. Results from Refs. [10] [42] other than fo B - Experi- 
mental quantities are taken from Ref . 1111 , except for f D + , 



det(stag + m) 

4 



1/4 



= det($ + m) 



(19) 



where the subscript denotes the number of flavors. The 
fourth-root can be built into chiral perturbation the- 
ory [ID]. In fact, in this context the j can be replaced 
by a free parameter, which is then fit. The fit yields 
0.28+ 0.03 [41], in excellent agreement with i. 

The two pri nciple methodological reasons why the er- 
ror in Eq. ([L2"]) is smaller than in Eq. J5]) is that Ref. [10] 
treats the charmed quark as a staggered quark [42] , us- 
ing a pseudoscalar density wi th an absolute normaliza- 
tion via a PCAC relation [13], and enabling an extrap- 
olation to the continuum limit. By contrast, Ref. [4] 
treats the charmed quark as a heavy quark [44] ; the cur- 
rent re qui res a matching factor computed in perturbative 
QCD [45], and the discretization effects are (conserva- 
tively) estimated with power-counting estimates |46] . 

To tackle charm on currently available lattices, the 
HPQCD Collaboration has developed a highly-improved 
staggered quark action (HISQ), first used to study char- 
momium [42] . Some of their other results are tabulated 
in Table [eT Especially noteworthy here is the value of 
fn+, which agrees with CLEO's later measurement [47] . 
Most effects that would bring Eq. Q12[) into better agree- 
ment with the measurements of fo s would also alter f D + , 
spoiling its agreement. 

For the it, K, and D + decay constants, both the chiral 
and continuum extrapolations are crucial. For the D s , 
however, the valence charmed and strange quarks ensure 
a mild chiral extrapolation. The continuum extrapola- 
tion turns out to be i nteresting: reading values for fu s 
off of plots in Ref. [10] , I have verified the continuum ex- 
trapolation and found that the slope in a 2 conforms with 
expectations of discretization effects of order a s a 2 m c A. 

The relevant portion of HPQCD's error budget is pre- 
sented in Table [2] Most of the row headings are self- 
explanatory, except for "scale r\ ," which is discussed be- 
low. The error budget is nearly complete, in my opinion, 



which is from Ref. 1471. Here A„ 



2m i 



d, s. 



Quantity 



Expt 



HPQCD units 



A, 
A s /A d 

u 
f K 

f D+ 



116.4+1.2 
758.7 
956.5 
1.261 + 0.002 
130.7+0.4 
159.8+0.5 
205.8+ 8.9 



111 + 5 
755 ± 14 
944+ 12 
1.252 + 0.015 
132 + 2 
157+2 
207 + 4 



MeV 
MeV 
MeV 

MeV 
MeV 
MeV 



Table 2. Error budget from Ref. 1101 . Entries in percent. 



Source 


Id. 


Sd 


fnjf, 


Statistics 


0.6 


0.7 


0.5 


Scale r\ 


1.0 


1.4 


0.4 


Continuum limit 


0.5 


0.6 


0.4 


Chiral limit 


0.3 


0.4 


0.2 


Adjust m s 


0.3 


0.3 


0.3 


Adjust m d © QED 


0.0 


0.1 


0.1 


Finite volume 


0.1 


0.3 


0.3 



more complete than many error budgets in the lattice- 
QCD literature. It does, however, fail to quote an uncer- 
tainty for quenching the charmed sea. This is surely a 
small effect, of order a s (A/m c ) 2 , but perhaps commen- 
surate with the \% errors included in Table |2~1 

I shall now discuss r\ in several steps, first motivating 
why it is used, then giving its definition and its value 
circa 2007. Being based on an expanding set of numeri- 
cal data, its value has now changed, so I discuss how it 
affects charmed-meson decay constants. 

Lattice gauge theory has a built-in ultraviolet cutoff — 
the lattice itself. The natural output is a dimension- 
less number, with physical dimensions balanced by pow- 
ers of the lattice spacing a. With a decay constant /, 
one computes af and then must introduce a definition 
for a. This is necessary not merely to quote a final result 
in MeV, but also to combine calculations at varying a, 
which is needed to understand the continuum limit. This 
is done by picking some fiducial mass M, and defining 
a = (aAf )i at /Mo X pt. This step eliminates one of the free 
parameters of QCD, namely, the bare coupling. 

To keep a long story short, no quantity is ideally suited 
to p lay the role of Af. A popular choice is 1/ri, defined 
via HH1HS] 



r 2 F(ri) 



(20) 



where F(r) is the force between two static sources of 
color, distance r apart. The advantages of r\ are that it 
is easy to compute in lattice QCD, and that it depends 
weakly on sea-quark masses and not at all on valence- 
quark masses. Then one can combine data from several 
lattices for nf = (ri/a)(af) in the chiral and continuum 
extrapolations. Other choices of M could complicate 
these steps. 

Of course, r± is unknown — it cannot be measured in 
the lab. It is inferred from the chiral and continuum 
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limit of other quantities. Reference [TU] used the value 

n = 0.321 ± 0.005 fm, (21) 

based on MILC's calculations of n/a [34] and HPQCD's 
own calculations of a(My( 2 s) — -^x(is)) [5Q]. The 1.6% 
uncertainty in r\ translates into a 1.0% uncertainty on 
fa, (c.f. Table I5T|) . because when r\ varies, the bare va- 
lence quark masses inside the D s do too. 

This rctuning when r\ changes has been studied by the 
Fermilab Lattice and MILC Collaborations (although 
not all details are as yet public). Since Ref. [1] was 
published, MILC has extended the ensembles to higher 
statistics, so Fermilab/MILC's decay constant analysis 
has continued, to reduce the total error. At Lattice 
2008 some of the discretization errors were brought un- 
der better control, leading to 249 ±11 MeV [5T , with 
(serendipitously) the same central value as Eq. ([5J ■ Ref- 
erences [1 EH] used n = 0.318 ± 0.007 fm [55] based on 
essentially the same input information as Eq. (|21[) . 

Meanwhile, however, evidence has begun to accumu- 
late that r\ should be smaller. Focusing on MILC's latest 
analysis of r\] v [53| . one has 



n = 0.3108 ± 0.0022 fm. 



(22) 



Rctuning the quark masses, this changes fo s to (prelim- 
inary, presented at Lattice 2009) 

f Ds = 260 ± 10 MeV Fermilab/MILC [55], (23) 

in which 4.2 MeV of the increase stems from the change 
in 7 *1, a nd the rest from other refinements of the analy- 
sis [52]. In other words, a shift down of 2.3% in ri has 
led to a shift up of 1.7% in fjj e . 

The HPQCD Collaboration has also incorporated the 
extensions of the MILC ensembles into its analysis of 
ri [54]. They find 



n = 0.3133 ± 0.0023 fm, 



(24) 



which is 2.4% lower than the value in Eq. (|2T]) . Although 
this suggests an increase in fu s of 3-5 MeV, one should 
keep in mind that HPQCD's calculations of fo s have 
proceeded to yet finer lattices. It seems prudent to wait 
for their own update, rather than applying a shift. 

Because both Fermilab/MILC and HPQCD use the 
same ensembles of lattice gauge fields, it is, or should 
be, a high priority to compute fo and fo 3 with other 
formulations of sea quarks. A promising development 
comes from the European Twisted-Mass Collaboration 
(ETMC), which has ensembles with rif = 2 over a range 
of sea-quark masses and lattice spacings (although not 
as extensive as MILC's with nj — 2 + 1). They find 
fo, = 244 ± 8 MeV [55], where the error stems from a 
thorough analysis of all uncertainties except the quench- 
ing of the strange quark. It is not easy to estimate this 
error reliably enough for averaging. Earlie r re sults with 
rbf = 2 obtained similar central values [56], or a bit 
higher [57] , albeit with larger error bars. 

It seems reasonable, then, to take as the current best 
estimate from lattice QCD, the weighted average of 
Eqs. (fT5 1) and P5 ]l : 



f Ds = 242.6 ± 2.9 MeV LQCD 2009, 



(25) 



which is 2.3cr lower than the average of measurements 
in Eq. (|18|) . The experimental statistical error continues 
to dominate this a, although if the central value of the 
lattice average were to increase by 3-5 MeV, the discrep- 
ancy would soften below 2a. 



4. New Physics 

The foregoing discussion makes clear that it is desir- 
able both for the experiments to improve further in pre- 
cision and for the latticc-QCD calculations to be con- 
firmed. Given the current status, it is conceivable that 
the tension will increase again to the point that it war- 
rants broad attention. With that in mind, this section 
provides some information on extensions of the Standard 
Model. 

The decays D s — > iv could be mediated by particles 
other than the Standard W, either by s-channel an- 
nihilation via another chargc-+l particle, by i-channcl 
exchange of a charge— 1-| particle, or by u-channel ex- 
change of a charge — | particle. All three kinds of par- 
ticle are popular enough in extensions of the Standard 
Model to have their own sections in the Review of Par- 
ticle Physics [HJ[5T]. The charge— 1-1 particle would be 
a W or a charged Higgs boson; the fractionally charged 
particles are known as leptoquarks. All would have a 
mass, presumably, at least as large as M\y- Their inter- 
actions can be parametrized by the effective Lagrangian 

C cS = V2G 1 a {sY15C){vl1^l) + V2G e P (s l5 c){9 L e R ) 
- V2G e v (s^c)(D L ^£ L ) + ^G%{- S c){v L i R ) 
+ V2G e T (sa^c)(D L a^£ R ), (26) 

where G l A and G l P appear in the leptonic-decay ampli- 
tude (j4]). The other interactions are likely to arise in 
non-Standard models, stemming from the chiral quan- 
tum numbers of the quarks and leptons. Nonzero G v 
would interfere with the leading Standard amplitude 
of the semileptonic decay D — > K£v 1 potentially mak- 
ing a significant change in the rate |58j . On the other 

hand, nonzero G s or G^ T would interfere with helicity- 
supprcssion, being visible only in an asymmetry of D — > 
K[iv after 10 7 or more events are recorded [55] . 

A W 1 alters the (semi)leptonic amplitude via Ga 
(Gy). Barring a carefully-built (i.e., finely-tuned) 
model, this is not a promising scenario [|5]. Many pop- 
ular charged Higgs models are also unpromising. Ref- 
erence [5] presents a charged Higgs model that could 
explain an excess of D s — > Iv events, but it predicts the 
same-sized excess in D + — > Iv. Now, however, this is dis- 
favored by the near-perfect agreement the most precise 
measurement of f D + g7j with lattice QCD g] \M [55J [55] . 

Leptoquarks, of several ilks, remain. Even here the 
charge-+| case is unpromising [9], owing to constraints 
from the lepton- flavor violating decays t — > fiss, where 
ss hadronizes to 4> or KK. This leaves the most promis- 
ing candidate to be an SU(2)-singlct, charged — | lep- 
toquark. This particle has the quantum numbers of a 
scalar down quark d, with an i?-violating interaction 



{n ce c L £ L - K qe V* s s L Vg L )d + n' cl c R t R d + H.i 



(27) 



where the superscript "c" denotes charge conjugation, 
and k and k' are coupling matrices. (With down squarks, 

the d field should take a family index, and n and n' yet 
another index.) Exchange of d generates Eq. (|56"|) with 



G 



G e v - K * ctKql V* s /AV2Ml 



r<i _ n£ — 

\JJ p UC 



V* s /4V2Mj = 2G{ 



(28) 
(29) 



Generalizations of Eq. (|27[) appear in non-Standard mod- 
els that arise in many contexts ESQ I6TJ ] |6"T ] 155 ] 153 ] [Ml [551 

M M M M EE ED EH E3 HfT^ 
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Fig. 1. Historical development of the fo g puzzle since June 
2005, in years since Ref. [4] was posted at arXiv.org. Cir- 
cles denote experimental measurements, from e + e — collisions 
at the T(45) in medium grey (orange); from e+e~ collisions 
at the Ds'^Dj*' threshold in dark grey (red); their weighted 
average is the medium grey (orange) dashed line with error 
band in light grey (yellow). Squares denote lattice-QCD cal- 
culations, with rij = 2 flavors in medium grey (cyan); with 
rif = 2 + 1 flavors in black (open symbols for conference re- 
ports); the 2 + 1 weighted average is the black dashed line 
with error band in grey. The discrepancy a is shown by the 
grey (green) solid line segments, labeled by the vertical axis 
on the right, (color online) 



The leptoquark possibility been examined in a broader 
context including constraints from dec ays of the D me- 
son, kaon, t lepton, and proton [75j . Reference [75j 
claims the se constraints make it difficult for the inter- 
action (|27|) to explain an enhancement in both B(D S — > 
\xv) and B(D S — ► tv) at the same time. 

5. Summary 

The developments of the fjj s puzzle are collected into 
Fig. [TTJ which presents the measurements and calcula- 
tions discussed above, and the time dependence of their 
respective averages. The time axis starts with the post- 
ing of Ref. [4], the first latticc-QCD calculation with 
2 + 1 flavors of sea quarks. The CLEO results of Ref. [5] 
(t = 1.8) and of Refs. [22l [23] (t = 3.5) each are averaged 
to reduce clutter. 

The discrepancy arose mostly with HPQCD's lattice- 
QCD calculation, and partly with the 2007 measure- 
ments from Belle and CLEO, rising to 3.8<7. In the past 
18 months, the tension in fo 3 has fallen to 2.3(7, from 

• CLEO 's new measurements of January 2009: — 0.8er; 

• Fermilab/MILC's (preliminary) update for Lattice 
2009: -0.13er; 

• HFAG's reinterpretation of BaBar's measurement: 
-0.67cr. 

• CLEO's new measurement of October 2009: +0.1a; 

The history of the discrepancy's significance is traced 
via the grey (green) piece-wise horizontal line and right 
vertical axis in Fig. [T7] Now. by the- way. if one follows 
Ref. [TT] and omits the BaBar result (as reinterpreted by 
HFAG), the discrepancy is 2.6a. 



The prospects for a resolution of the /rj s puzzle arc 
good — whether the tension goes away completely or 
tightens again. BaBar is measuring the absolute branch- 
ing ratio, and Belle plans to update its analysis with 
higher statistics. In a few years, BES 3 will measure 
D s — > iv in threshold production, similarly to CLEO, 
with a target uncertainty of 1% [76]. Several lattice- 
QCD collaborations now have enough n/ = 2 + 1 ensem- 
bles to carry out a useful calculation of fo B ■ The MILC 
Collaboration has started to generate ensembles with 
n,/ = 2 + 1 + 1 sea quarks with the HISQ action, where 
the fourth sea quark is charm. Similarly, the ETMC has 
embarked on a project with n/ = 2 + 1 + 1 twisted-mass 
Wilson sea quarks. Even if the puzzle dissipates, D and 
D s leptonic decays will be useful for constraining exten- 
sions of the Standard Model [75] , 
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